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Weak-lensing distortions of the cosmic-microwave-background (CMB) temperature and polariza¬ 
tion patterns can reveal important clues to the intervening large-scale structure. The effect of leasing 
is to deflect the primary temperature and polarization signal to slightly different locations on the 
sky. Deflections due to density fluctuations, gradient-type for the gradient of the projected gravita¬ 
tional potential, give a direct measure of the mass distribution. Curl-type deflections can be induced 
by, for example, a primordial background of gravitational waves from inflation or by second-order 
effects related to leasing by density perturbations. Whereas gradient-type deflections are expected 
to dominate, we show that curl-type deflections can provide a useful test of systematics and serve 
to indicate the presence of confusing secondary and foreground non-Gaussian signals. 
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I. INTRODUCTION 

The theory of weak gravitational leasing (“cosmic 
shear”) of the cosmic microwave background (CMB) by 
large-scale mass inhomc^eneities has now been studied 
extensively niimEn The distortions produce 
a unique non-Gaussian signal in the CMB temperature 
and polarization patterns. The distinction with the pri¬ 
mordial Gaussian pattern can be exploited to determine 
various properties of the leasing sources. Since inhomo¬ 
geneities in the intervening mass field are expected to 
provide the dominant source of leasing, likelihood meth¬ 
ods and quadratic estimators have been developed to 
map the projected mass distribution. Such reconstruc¬ 
tion will be important not only to study the matter power 
spectrum, but also to reconstruct the primordial polar¬ 
ization signal from inflationary gravitational waves 0 . 

To linear order in the density-perturbation amplitude, 
leasing by mass fluctuations results in a deflection an¬ 
gle that can be written as the gradient of a projected 
gravitational potential (i.e., the deflection angle is a lon¬ 
gitudinal vector held in the plane of the sky). Algorithms 
to reconstruct the mass distribution hence measure only 
this longitudinal component of the deflection angle. A 
deflection angle that can be written as a curl—a gradient- 
free or transverse-vector held—can be produced through 
lensing by gravitational waves or through lensing by mass 
fluctuations to second order in the density-perturbation 
amplitude. Since the curl-type deflection is expected to 
be significantly weaker than the gradient-type (as dis¬ 
cussed further below), measurement of the curl power 
spectrum has been used as a test of systematic artifacts 
in the cosmic-shear maps that have been produced with 
measurements of shape distortions to high-redshift galax¬ 
ies For galaxy cosmic-shear maps, this curl com¬ 

ponent is measured by simply rotating each galaxy image 
by 45° [ISini. 

For cosmic-shear of the CMB, one cannot simply ro¬ 
tate the temperature pattern at each point on the sky. 


However, there is indeed a method to reconstruct the curl 
component of the deflection angle that is directly analo¬ 
gous to that for reconstructing the gradient component 
[ 3 . In this paper, we show that the cosmological curl 
signals are expected to be small, and that measurement 
of the cosmic-shear curl component can thus be used as 
a diagnostic for systematic artifacts, unsubtracted fore¬ 
grounds, and/or primordial non-Gaussianity. 

This paper is organized as follows. In Section m we 
introduce our formalism and present expressions for the 
weak-lensing corrections to the CMB power spectra from 
both gradient and curl modes of the deflection field. In 
Section ITTTI we discuss quadratic estimators for the deflec¬ 
tion field for both the gradient and curl components, and 
we discuss the orthogonality of these estimators or Alters. 
Section nm discusses cosmological sources for a curl 
component, first gravitational waves and then second- 
order density perturbations. Section El presents results 
of our calculations. Section EH discusses the use of a 
curl reconstruction as a diagnostic for primordial non- 
Gaussianity, unsubtracted foregrounds, or systematic ar¬ 
tifacts. A few concluding remarks about the applications 
of our work are given in Section IVIII 

II. EFFECT ON CMB POWER SPECTRA 

The effect of weak lensing on CMB anisotropies is a 
non-linear remapping of temperature and polarization 
fluctuations. In the case of temperature anisotropies on 
the sphere, this remapping can be expressed as 

0(n) = 0[n -I- V(/(n) -|- V x H(n)], (1) 

where 0 is the observed temperature pattern, n = (0, (p) 
is the position on the sky, and 0 is the original unlensed 
pattern. Here, the gradient Vc/ has components dip in 
the plane of the sky, and the “curl” V x O has components 
eijdjUl, where Cij is the antisymmetric tensor. The map¬ 
ping involves the angular gradient of the projected gravi- 
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tational potential (j) due to density perturbations, and the 
curl of some other function fl, to be discussed further be¬ 
low. The convergence n and the image rotation to that 
usually arise in gravitational leasing of discrete sources 
are given by /tfn) = —^V^(j){n) and Ct’(h) = —iV^n(n), 
respectively [Ij, [l3 ■ In fh® limit of weak deflection, the 
remapping in Eq. m can be expressed in Fourier space 
as 

0 ( 1 ) = J dn 0 (n)e-*' " 

= e{i)- 1 (2) 

under the flat-sky approximation, where 

L^(l,l') = fl(l-1') [(1-1') X !']• (4) 

Note that the curl component has a two-dimensional 
cross-product, a ninety-degree rotation between Fourier 
components which we denote as x, following Ref. 0 , 
whereas the gradient has a dot-product. Strictly speak¬ 
ing, there is an additional term quadratic in (j) and fl, 
respectively, that is required to obtain the lowest-order 
cosmic-shear corrections to the power spectrum; see, e.g., 
Eq. (3) in Ref. Q, and for simply replace the dot 
products therein by cross products. For economy, we do 
not reproduce those expressions here, but they are in¬ 
cluded in our numerical work. 

I 

cr = ii-Ricr+i/|i5qr4,[(i- 
dr = ii-/iicr+i/0qf4,i(i- 


In the above, the operator 0 is a dot product when X = (j) 
and a cross product when X = Q. Note that inclusion 
of the R correction, which was neglected in Ref. |^, is 
required to obtain a correction that is complete to lowest 
nonvanishing order in or fl. We will present numerical 
results for these power spectra in Section El after dis¬ 
cussing the power spectra and in Section HV Al 


III. RECONSTRUCTION OF THE 
DEFLECTION FIELD 

So far we have discussed the corrections to the tem¬ 
perature/polarization power spectra due to weak gravita¬ 
tional leasing. However, what is perhaps more interesting 
is that lensing induces characteristic non-Gaussian signa¬ 
tures in the temperature/polarization pattern. Measure¬ 


The temperature-anisotropy power spectrum is 

cf = [1 - i?] cf + J |^qtp|Qf^[(i - ii) © ii]". 

(5) 

Here, is the power spectrum of either lensing po¬ 

tentials related to density fluctuations or the rotational 
component, while i? is a multiplicative correction 0{4>^) 
that can be obtained from Eq. (8) in Ref. 0 (again re¬ 
placing a dot product by a cross product for H). 

In addition to temperature anisotropies, lensing also 
modifies the polarization. We follow the notation in 
Ref. 0 , and then the remapping of the polarization un¬ 
der lensing is 

±X(h) = ±X[n -I- V0(n) -|- V x H(n)] (6) 

where ^X = Q ± iU. Since the Stokes parameters are 
not rotationally invariant, we write them in terms of the 
rotational invariants E and B 0 which are defined in 
Fourier space through ±X(^) = [i?(/)=F*S(0]e^^*‘^, where 
ip is the phase angle of 1. Then, the observed polarization 
is 

±X{1) = ±X{1) - J . 

_ ^ (7) 

Following Ref. |2| , the lensed polarization power spectra 
can now be expressed in terms of for X = {(j), H) 

and the unlensed CMB spectra as 

li) 0 cos 2(2^,J + sin2(2^,J], (8) 

li) © + cos2(2^q)]. (9) 


ment of these non-Gaussianities can be used to map the 
deflection-angle as a function of position on the sky, and 
thus to infer the projected potentials (j) and H. 

We now extend the quadratic estimators that have 
been proposed to reconstruct the gradient component of 
the deflection field dillli to the case of a curl com¬ 
ponent @. In Fourier space, the quadratic estimator can 
be written for X = (/> or X = H as 

^(1) = / li)0(li)0(l - li), (10) 

where is a filter that acts on the GMB tempera¬ 
ture field subject to the demands that (< 7 ^( 1 )) = (()(1) and 
(H(l)) = H(l). The filters that optimize the signal-to- 
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noise are 


W^{lAi)=Ni 


X 






( 11 ) 


f dH, {(l0li)Ce + [l0(l-lO]C,t,,|P 

"7 ( 2^)2 

( 12 ) 

Here, Cf is the total temperature power spectrum and 
can be written as a sum of the lensed power spectrum, 
foregrounds, and detector noise: Cf = 

Filters similar to these can be written down for the other 
quadratic combinations of polarization and temperature. 
We do not write them out explicitly here as they can be 
derived easily from published expressions in the litera¬ 
ture. 

In practice, one determines each Fourier mode ^(1) 
[or H(l)] by taking an appropriately weighted average of 
all combinations of temperature (or polarization) Fourier 
modes li and I 2 that sum to li -b I 2 = 1. The only differ¬ 
ence between the reconstruction of gradient versus curl 


J 


modes is whether to weight these combinations by a dot 
product li • I 2 or by a curl li x I 2 . The quantity is the 
noise, the variance with which each Fourier mode (/)(!) or 
H(l) can be reconstructed. Thus, when the (j) or Q power 
spectrum is measured with these quadratic estimators, 
the power spectra of the estimators will be, 

= (27r)2^(l + . (13) 

Of course, if the power spectra Cf and C® are known, 
then the noise can be calculated independently and sub¬ 
tracted to yield the desired 4> or O power spectra. 

The orthogonality of the weightings in the filters for O 
and (j) suggests that if we have a deflection field that is a 
pure gradient, then the application of the curl filter will 
give zero, and vice versa. Although this is approximately 
correct, it is not precisely true, as we now show. Con¬ 
sider a deflection field that is a pure gradient; i.e., it is 
described in terms of nonzero ^(1), with H = 0. Suppose 
now that we measure ft: taking Eq. cnj with X = ft, 
the only possible source in the temperature field is due 
to the gradient, whereby 


(H(l)) = 


<^(1) / ^VF^(l,li){(Mi)C® + [l-(l-li)]C|®_i^l} 

dH, {(lxli)C® + [lx(l-b)]C|ti,|} 




( 2^)2 


0 /^t 

"'-X 7 11 1 I 


{(Mi)C® + [l-(l-li)]C'®_q|} . (14) 


Despite the fact that the filter is designed to select out 
only the curl contribution, close inspection of this integral 
shows that it is not precisely zero. Note, however, that if 
Cf is a pure power-law — thus, also Cf is a power-law 
— then the integral would vanish identically. In other 
words, the (p and ft filters are orthogonal only to the 
extent that Ci behaves as ^ The departure from or¬ 
thogonality is due to the presence of bumps and wiggles 
in the CMB anisotropy power spectrum. The departure 
from a power-law spectrum also prevents the construc¬ 
tion of precisely orthogonal filters to separate the two 
modes exactly. 

The integrand in Eq. m is nonzero only for values 
of h at which C® and C'®_j^| departs from the power- 
law. The departure is not signihcant except when li or 
1 —111 enter the damping tail of the anisotropy spectrum. 
This contribution, however, is suppressed by finite angu¬ 
lar resolution of CMB experiments. We therefore expect 
that the integral in Eq. da should be small, even if it is 
not precisely zero. A numerical evaluation conhrms this 
argument; we have found that the expression evaluates 
to well below I0“^° (p{l)N[^ for I values up to 5000. We 
therefore conclude that the reconstruction can be consid¬ 


ered to be effectively orthogonal. 

This leads us to another point. Cosmic shear, either 
through cj) 01 ft, leads to a correction to the observed 
CMB temperature that can be written as (V6) • (Vt/f) 
or (V0) • (V X ft). Suppose, however, that some other 
process lead to a correction of the form 0(n)/(n), where 
/(n) is some function of position on the sky. For ex¬ 
ample, consider the Sunyaev-Zel’dovich effect m The 
thermal effect can be subtracted to a large extent through 
multifrequency observations. However, the kinetic-SZ ef¬ 
fect has the same frequency dependence as primordial 
fluctuations and is therefore indistinguishable. Suppose 
that there is thus some unsubtracted SZ contribution 
to the measured CMB fluctuation. Then this will pro¬ 
vide an angle-dependent multiplicative correction to the 
primordial temperature. Something similar (though not 
precisely so) may occur through exotic phenomena such 
as primordial non-Gaussianity or a spatially-varying 
fine-structure constant Ea, for example. On the other 
hand, non-uniformities in the instrumental gain may also 
mimic such an effect. A quadratic estimator can be con¬ 
structed for /(n), simply by removing the 1 vector depen¬ 
dences in Eqs. CD and CD- Again, the estimator for / 
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will be close to orthogonal to those for (p and fl. However, 
the orthogonality will not be precise, and if there is a sig- 
nihcant /(n), then it will show up in a reconstruction of 
(j) and to a similar level in fl. 


IV. COSMOLOGICAL CURL SOURCES 

A. Primordial Gravitational Waves 

Our first example of a curl deflection is a background of 
gravitational waves from inflation. Refs. [i3,ll3 showed 
that gravitational waves can act as gravitational lenses, 
and Refs. MM showed that leasing by gravitational 
waves gives rise to a curl component in the deflection an¬ 
gle. Suppose there is a gravitational wave with amplitude 
hij (more precisely, the transverse traceless tensor part 
of the metric perturbation), and suppose further that 
we choose our line of sight to be (near the) z direction. 
Then, H oc tudkhzi- For example, if the gravitational 
wave propagates in the y direction, then H oc dyhzx and 
the deflection is in the 9^ direction with 59^ oc hzx- 
Of course, the total deflection is an integral of all the 
deflections along a line of sight, and for arbitrary line of 
sight, the rotation is [l^ . 

w(n) = —n • [Vfi X r(n,rs)] 

2 Ts 

= - i dr'[n.(VxH).n] , (15) 

where H is the transverse (V • H = 0), traceless (Tr H = 
0 ), tensor metric perturbation representing gravitational 
waves. In the above, r is the radial distance from 
the observer and the surface of last scattering is at rg- 
The gravitational-wave amplitude obeys a wave equation 
which, ignoring the presence of anisotropic stress from 
neutrinos and other relativistic species at early times, 
takes the form H — V^H -|- 2(d/a)H = 0, where the dot 
denotes derivative with respect to conformal time. We 
express the solution to this equation in the form of a 
transfer function, r(T)(k;, r), whereby the Fourier ampli¬ 
tude evolves as H (k, r) = H (k)T(x) {k, r) and H (k) is the 
initial amplitude of the wave. In a purely dust-dominated 
universe (appropriate for the long-wavelength modes of 
relevance here, which come into the horizon at late time), 
T(T)(fc,r) = 3ji(Ar)/(fcr) (Ref. [H presents more precise 
expressions, but they are not relevant for the calculation 
here). Assuming isotropy, the three-dimensional spatial 
power spectrum of initial metric fluctuations related to a 
stochastic background of gravitational waves is 

(%(k)id(b)(k')) = (27r)3p(x)(fc)<5„ d(3)(k-k0 , (16) 

where the two linear-polarization states of the gravita¬ 
tional wave are denoted by {i,j)- In standard inflationary 
models, the primordial fluctuation spectrum is predicted 


to be 


P(^T){k) = (17) 

Inflationary models generally predict that tit ~ 0 while 
the ratio of tensor-to-scalar amplitudes, r = At/A s, is 
now constrained to be below 0.36 M We will use the 
upper limit allowed when calculating the inflationary- 
gravitational-wave (IGW) contribution. 

Taking the spherical-harmonic moments of Eq. 113, 
the angular power spectrum of the rotational component 
is 




1 


21 + 1 




m— — l 


k^dkP(T){k)\T/^{k,rs)\\ (18) 


where 


rr(fc,r) 



(19) 


For comparison, the gradient components of the deflec¬ 
tion angle involve the projected density perturbations 
along the line of sight to the last-scattering surface, 

(^(m) = —2 / dr^ —^<i)(r, thr), (20) 

Jo [rrs] 

where $ is the potential associated with the large-scale 
mass distribution. The angular power spectrum of these 
projected potentials are 




I 


21 + 1 
2 
TT 




k^dkP^s)ik) Tf{k,rs) 


( 21 ) 


where P(^s){k) is the power spectrum of density perturba¬ 
tions, including the transfer function, and 


T^ik,rg) = -3a 



,Gir')rs -r' 

dr -j-A— - 3i{kr ) 

a[r') r+s 


( 22 ) 

with the growth of matter fluctuations given by G(r), 
and a{r) is the scale factor. Here, we ignore the met¬ 
ric shear at the surface of last scattering by the same 
background of waves. As discussed in Ref. [i3|, the curl 
component due to intervening deflections from the grav¬ 
itational background is small compared to the gradient 
component and the inclusion of metric shear only leads 
to a further cancellation. Hence, our forthcoming pro¬ 
posal that the curl component should be considered as a 
test of systematics will not be affected adversely. 
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FIG. 1: Lensing-deflection power spectra. Here, we show the 
gradient component from density perturbations (top curve), 
the curl component from inflationary gravitational waves 
(dashed curve labeled ‘IGWs’), and the curl component from 
second-order density perturbations, (dot-dashed curve). We 
have taken the maximum IGW amplitude consistent with the 
current upper limit to tensor-to-scalar ratio IT^ . 


B. Second-order density perturbations 

Gravitational leasing by density perturbations can give 
rise to a curl component once we go to second order 
in the projected potential (j). To see this, we first re¬ 
view the lowest-order effect. Suppose there is a lens 
at a distance ri along the line of sight. The deflec¬ 
tion by this lens is S0i oc 9i$i, where is the grav¬ 
itational potential (not the projected potential) at ri. 
The lowest-order deflection will therefore be written as 
the sum of gradients perpendicular to the line of sight. 
To second order in <&, there can be deflection by two 
lenses at different distances, ri and r2, along the line 
of sight. The deflection by the first lens the ray en¬ 
counters is oc 9j$i, and the deflection after encounter¬ 
ing the second lens is oc (9i9j$2)(9j$i) (this follows 
from the discussion, e.g., in Section 3 of Ref. 53)- To 
see that this has nonvanishing curl, we take the curl: 
eikdk[{didj^ 2 )idj^i)] oc eik{didj^ 2 ){djdk^ 2 ) which does 
not generally vanish. 

The full calculation of the curl power spectrum is then 
lengthy but straightforward, and it is discussed in the 
context of galaxy-based weak leasing surveys in Ref. 0 
They are explored in the context of CMB leasing in 
Ref. 13. We do not repeat the derivation but refer the 
reader to Ref. 0 for details. 


V. CALCULATIONS AND RESULTS 

A comparison of the gravitational-wave and density- 
perturbation (to second order) curl signals to the gra¬ 


dient lensing signal is shown in Fig. ^ [Note that the 
anisotropy spectra for the lensing convergence and ro¬ 
tation are related to the gradient and curl by = 
{l^/4)Cf^ and Here, we push the 

gravitational-wave amplitude to the maximum allowed 
by current data. A useful measure of the relative impor¬ 
tance of the two components is the rms deflection angle 
on the sky given by 9"^^^ = j[(P\/. In the 
case of density perturbations (for X = cj)), 9rnis = 7x 10“^ 
or roughly 2.5 arcmins. The angular coherence scale, 
where the rms drops to half its peak value, is about a 
degree. In the case of the strongest gravitational wave 
background, the deflection angle is 0rms = 7 x 10“® or 
0.25 arcmins, but the angular coherence scale is a few 
tens of degrees. (Note that the y axis in Fig. ^is l^Ci, so 
the power spectra plotted there are in fact very rapidly 
falling with 1. The rms deflection angle is thus fixed pri¬ 
marily by the low I’s.) In the case of second-order curl 
corrections, the coherence scale is similar to that of den¬ 
sity perturbations but the amplitude is smaller by at least 
four to five orders of magnitude. The resulting correc¬ 
tions to CMB anisotropies trace that of the density field, 
but with a similar reduction in the overall amplitude. 

The effects on the CMB anisotropy spectra are sum¬ 
marized in Fig. 121 In the case of temperature, the 
gravitational-wave-lensing correction is at least two or¬ 
ders of magnitude below the temperature fluctuations 
associated with the angular displacement corrections due 
to the density field. We also summarize our results for 
the case involving polarization anisotropies. In accord 
with Ref. 0 , we conclude that the correction resulting 
from the curl component is negligibly small. 

We now turn to the reconstruction of the cosmic-shear 
pattern with quadratic estimators for cj) and H. Fig. |21 
shows the errors in the reconstruction for a hypothetical 
CMB experiment with a resolution of an arcminute and 
a noise-equivalent temperature of 1 yK sec^/^ over one 
year of integration. We show the reconstruction for both 
temperature maps (top lines), and for the EB quadratic 
combination which was shown in Ref. 0 to be the best 
combination to extract lensing information from CMB 
data. Note that one generally reconstructs the gradient 
and curl components of the deflection field with roughly 
the same signal-to-noise ratio. The gradient component, 
however, dominates since it is sourced by the large-scale 
mass distribution, while the curl component is subdomi¬ 
nant given that the amplitude of the tensor contribution 
to the CMB quadrupole is limited by current CMB data 
to be less than 30% of that due to scalar perturbations. 


VI. THE CURL AS A DIAGNOSTIC 

Is it useful to reconstruct the curl component as there 
is virtually no signal? Here, we suggest that a recon¬ 
struction may be useful to identify non-Gaussian sig¬ 
nals, both due t o p rimary effects, such as primordial 
non-Gaussianity |2nl| or perhaps variable fine-structure 
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FIG. 2: The lensing modification to CMB power spectra for density perturbations and for gravitational waves. Left: Tem¬ 
perature fluctuations. The top curve is the primordial power spectrum. The middle curve is the (additive) contribution to 
the temperature power spectrum from lensing by density perturbations, and the lower curve is for lensing by gravitational 
waves, assuming the maximum IGW background consistent with current GMB bounds. Note that for the lower two curves, it 
is \Ci — Ci\ which is plotted. The negative contribution of the R term in Eq. 0 allows for Ci — Ci to become negative. Since 
the coherence scale for is so small, the lensed power spectrum reflects closely at high I the primordial temperature power 
spectrum. Right: Polarization. These curves are: (1) the top (solid) curve is the primordial E-mode power spectrum; (2) the 
next-highest (solid) curve is the lensing correction to the primordial E-mode power spectrum by density perturbations; (3) the 
top dashed curve is the B-mode power spectrum resulting from cosmic-shear conversion of E-modes by density perturbations; 
(4) the lower dashed curve is the B-mode power spectrum from lensing by foreground IGWs; and (5) the lowest solid curve 
is the E-mode power spectrum resulting from lensing by foreground IGWs. The dot-dash curve is the primordial B power 
spectrum from the maximal IGW background allowed by current constraints. 



FIG. 3: Reconstructed unbinned noise spectra using the 
temperature-temperature quadratic estimator and the EB po¬ 
larization combination. We show noise for both grad and curl 
modes. For reference, we also plot the power spectrum of the 
deflection angle corresponding to the gradient component. 


constant and secondary anisotropies. As discussed 
above, although the thermal SZ effect can can be removed 
largely from multifrequency data [IJ, the kinetic-SZ ef¬ 
fect cannot. Such an unsubtracted secondary contribu¬ 
tion to CMB fluctuations will give rise to an additional 


noise-bias term in Eq. whereby 

(A(l)A(l')) = {2^f5{\ + l')(Cf^ + S^). (23) 

When secondary non-Gaussianities are not properly ac¬ 
counted for, the additional noise-bias term takes the form 

sf = {Nn" 

f dH, f dH, {(l©li)C® + [l©(l-li)]ce_,^|} 

^ J {2nyJ (27r)2 

{(l©l2)Ce + [l©(l-l2)]C®_,^l} 

0/^t /^t 

X (0*(li)0*(l-li)0*(l2)0^(-l-l2)), (24) 

where (0*(li)0"(l - li)0*(l2)0"(-l - b)) is the four- 
point correlator of the contaminant foreground or pri¬ 
mordial non-Gaussianity with its anisotropy written in 
Fourier-space as 0®(1). This correlator can be decom¬ 
posed as 

(0*(ll)0^(l-ll)0^(l2)0*(-l-l2)) 

= 2Q7C|r_,p<5(li + 12 ) 

+r^(li,l-li, 12,-1-12). (25) 

The Gaussian piece leads to a noise bias 

= {Nn" 
























7 



^ 21 ^ {(l0b)Ce+ [!©(!- 


(27r)2 

nqL-,.i 

2 


X 


2Q7C' 


ss 


(26) 


which can be absorbed into with a proper definition 
of the normalization factor, and where also include 
foregrounds and secondary power spectra. However, the 
non-Gaussian nature of the foreground cannot be ignored 
and this results in a bias that cannot be removed by a 
renormalization. This noise is, 




(2nr 

{(l©li)ce + [l©(l-li)]C|t,,|} 
{(l©l2)Ce + [l©(l-l2)]C|t,,|} 

X T«(li,l-li, 12,-1-12). 


(27) 


The angular dependence of the foreground trispectrum 
is important: if the trispectrum were to depend on the 
length of the vectors alone, the averaging would result in 
significant suppression of this noise bias. 

In the presence of an additional secondary trispectrum, 
for lensing reconstruction, 

(<^(1)<^(1')) = (27r)2<5(l + 10(67/-^ + Nf + Sf ). (28) 

While Nf can be established based on noise properties, 
one cannot separate the signal of interest from the 
confusion Sf. Such a situation has already been ob¬ 
served, for example, in numerical simulations of the CMB 
lensing reconstruction process where in the presence of a 
kinetic-SZ component an additional noise bias was sug¬ 
gested . The presence of such a noise bias was readily 
detectable given that the input mass spectrum, or alter¬ 
natively the cosmology, was known a priori. In the real 
case, one is interested in measuring these quantities from 
the mass spectrum determined from CMB lensing. Thus, 
the presence of a noise bias cannot easily be established 
since the bias is degenerate with other unknowns. 

The curl component, however, provides a useful 
method to establish the presence of such a noise bias 
which can be used to correct the gradient component or 
to allow for an accounting of the bias when cosmological 
parameters are measured. This follows from the fact that 
all signals of interest in the curl component are negligibly 
small such that the resulting reconstruction only leads to 

(H(l)H(l')) = {27TfS{\ + l'){Np + Sf^). (29) 

Since TVp, like Nf can be predicted from the measured 
CMB power spectra, any excess noise in the reconstruc¬ 
tion will suggest either systematics or the presence of an 


additional non-Gaussian signal that is contributing via 
5'p. Even though the excess noises in the gradient and 
curl of the deflection field are different—i.e., S'p vs Sf — 
the origin of the excess noise could very well be the same 
with the only differences resulting from variations in the 
two filters for the two modes. In general, any detection 
of excess noise in the curl component should suggest a 
bias in the gradient component. Since filter shapes are 
known a priori^ one should be able to establish some 
estimate on the expected excess noise in the gradient, 
given the excess noise in the curl. If this excess noise 
is significant, then the dominance of a systematic effect 
in the reconstruction is clearly established. Currently, 
there is no mechanism to either estimate or establish the 
presence of a systematic noise component in the CMB 
lensing analysis. Thus, we suggest that the curl compo¬ 
nent be used as a monitor of systematic effects and to 
understand if the Cf^ reconstruction is affected through 
Si by non-Gaussian secondary effects and foregrounds. 

In general, we do not expect effects such as primordial 
non-Gaussianity to be a significant concern for lens¬ 
ing reconstruction of the deflection-potential statistics. 
Given the noise levels to the reconstruction, as shown 
in Fig. 0 one can establish the minimum amplitude for 
which systematic effects or additional noise biases, as de¬ 
scribed by Si , can be detected via 



Here ai = y^2/(2l -|- l)Ni, under the assumption of no 
signal in the curl component. Using the estimated noise 
levels, from the EB combination of polarization maps to 
reconstruct the curl component, for example, one can es¬ 
tablish systematic effects down to a level of 0.1% from 
the amplitude of the potential-fluctuation power spec¬ 
trum. In the cosmic-shear simulations of Ref. noise 
biases at the level of 30% or more were found. We surmise 
that some of this may be due to conversion of kinetic-SZ 
corrections to a deflection angle, although this probably 
does not account for all the excess noise. A study of the 
curl component may help clarify the nature of such noise 
biases in the simulation. 


VII. SUMMARY 

Lensing by gravitational waves is expected to give rise 
to a cosmic-shear pattern where the deflection angle has 
a curl (or transverse-vector) component, as opposed to 
the curl-free pattern expected by cosmic shear by density 
perturbations (to linear order in the perturbation ampli¬ 
tude). To second order in the perturbation amplitude, 
lensing can also give rise to a curl component. For a pri¬ 
mordial background of gravitational waves from inflation 
with a normalization given by the current upper limit 
to the tensor-to-scalar ratio, the corrections to the CMB 










power spectra are generally two orders of magnitude be¬ 
low those of the cosmic shear due to density perturba¬ 
tions, and the curl component from higher-order lensing 
effects is also small. The curl component can be recon¬ 
structed with quadratic estimators analogous to those 
developed to measure the gradient component. Given 
the small signal expected in the curl, this component can 
potentially be used as a probe of systematic effects and 
foregrounds for next-generation CMB experiments 
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